Stable axisymmetric SPH formulation with no axis singularity          下载免费PDF全文

Jian Wang  Dave Chan  Jiahe Zhang 《国际地质力学数值与分析法杂志》2016,40(7):987-1006

The axisymmetric formulation of the governing equations for geomechanics in the framework of smoothed particle hydrodynamics (SPH) is presented in this study. Two forms of SPH discretization for the motion equations, which are labeled as form I and form II, are proposed, and the methods to compute the hoop stress and strain terms including hoop strain rate and the acceleration introduced by the hoop stress are compared. To avoid possible singularity problem near the axis of symmetry, a perfectly smooth contact along with ghost particles are applied to prevent the real particles from overly approaching the axis of symmetry to remove this potential singularity. In addition, the Mohr–Coulomb constitutive model is implemented into the SPH formulation in describing soil behavior. Four numerical tests are carried out to validate and compare the accuracy and stability of the proposed algorithms, and their results are compared with analytical solutions and results from FEM analysis. The performance in these comparisons suggests that SPH II with hoop terms computed through direct hoop method is more stable than the others, and the adoption of contact for the symmetric axis is efficient in eliminating the singularity problem. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

A nonstationary extreme value distribution for analysing the cessation of karst spring discharge          下载免费PDF全文

Yan Liu  Yonghong Hao  Yonghui Fan  Tongke Wang  Xueli Huo  Youcun Liu  Tian‐Chyi J. Yeh 《水文研究》2014,28(20):5251-5258

The effects of climate change and population growth in recent decades are leading us to consider their combined and potentially extreme consequences, particularly regarding hydrological processes, which can be modeled using a generalized extreme value (GEV) distribution. Most of the GEV models were based on a stationary assumption for hydrological processes, in contrast to the nonstationary reality due to climate change and human activities. In this paper, we present the nonstationary generalized extreme value (NSGEV) distribution and use it to investigate the risk of Niangziguan Springs discharge decreasing to zero. Rather than assuming the location, scale, and shape parameters to be constant as one might do for a stationary GEV distribution analysis, the NSGEV approach can reflect the dynamic processes by defining the GEV parameters as functions of time. Because most of the GEV model is designed to evaluate maxima (e.g. flooding, represented by positive numbers), and spring discharge cessation is a ?minima’, we deduced an NSGEV model for minima by applying opposite numbers, i.e. negative instead of positive numbers. The results of the model application to Niangziguan Springs showed that the probability of zero discharge at Niangziguan Springs will be 1/80 in 2025, and 1/10 in 2030. After 2025, the rate of decrease in spring discharge will accelerate, and the probability that Niangziguan Springs will cease flowing will dramatically increase. The NSGEV model is a robust method for analysing karst spring discharge. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

Temporal variation of soil erosion resistance on sloping farmland during the growth stages of maize (Zea mays L.)     

Rui Ma  Zicheng Zheng  Tingxuan Li  Shuqin He  Xizhou Zhang  Yongdong Wang  Huagang Huang  Daihua Ye 《水文研究》2021,35(9):e14353

Maize growth has great effects on soil properties and thus likely induces the changes in soil erosion resistance on sloping farmland. However, temporal variation of soil erosion resistance during the growth stages of maize is still unclear in the mountainous yellow soil area where maize is the dominant crop. In this study, four maize plots (MP) and four bare land plots (CK) were conducted to investigate soil erosion resistance, and multiple indicators of soil erosion resistance were measured including the total soil anti-scourability (TAS), mean weight diameter (MWD), soil erodibility K factor and soil shear strength (SH). A comprehensive soil erosion resistance index (CSERI) was employed to quantify the temporal variation of soil erosion resistance during the growth stages of maize (seedling stage, SS; jointing stage, JS; tasselling stage, TS; maturing stage, MS). The results showed that TAS, MWD, SH increased significantly with maize growth and SH decreased when at MS. But K factor decreased significantly over time. CSERI increased significantly during the growth stages of maize and the CSERI of JS, TS, MS increased on average by 74.72, 180.68 and 234.57% than that of SS. Compared to CK, CSERI of MP increased by 49.90, 66.82, 55.60 and 38.61% during the growth stages of maize. The temporal variation of soil erosion resistance was closely related to the changes in maize cover, maize roots and soil organic carbon. The findings demonstrated that it is necessary to consider the temporal variation of soil erosion resistance in the mountainous yellow soil area.  相似文献   

Impacts of groundwater flow on the evolution of a thermokarst lake in the permafrost–dominated region on the Qinghai–Tibet plateau     

Jinlong Li  Wei Wang  Xianmin Ke  Jianhua Li  Yating Yu  Zeyong Gao  Fujun Niu 《水文研究》2021,35(12):e14443

As a result of global warming induced permafrost degradation in recent decades, thermokarst lakes in the Qinghai–Tibet plateau (QTP) have been regulating local hydrological and ecological processes. Simulations with coupled moisture–heat numerical models in the Beiluhe basin (located in the hinterland of permafrost regions on the QTP) have provided insights into the interaction between groundwater flow and the freeze–thaw process. A total of 30 modified SUTRA scenarios were established to examine the effects of hydrodynamic forces, permeability, and climate on thermokarst lakes. The results indicate that the hydrodynamic condition variables regulate the permafrost degradation around the lakes. In case groundwater recharges to the lake, a low–temperature groundwater flow stimulates the expansion of the surrounding thawing regions through thermal convection. The thawing rate of the permafrost underlying the lake intensifies when groundwater is discharged from the lake. Under different permeability conditions, spatiotemporal variations in the active layer thickness significantly influence the occurrence of an open talik at the lake bottom. A warmer and wetter climate will inevitably lead to a sharp decrease in the upper limit of the surrounding permafrost, with a continual decrease in the duration of open talik events. Overall, our results underscore that comprehensive consideration of the relevant hydrologic processes is critical for improving the understanding of environmental and ecological changes in cold environments.  相似文献   

Chemical weathering and CO2 consumption in the glaciated Karuxung River catchment,Tibetan Plateau     

Fan Zhang  Xiong Xiao  Lijie Wang  Chen Zeng  Zhengliang Yu  Guanxing Wang  Xiaonan Shi 《水文研究》2021,35(8):e14330

Climate factors play critical roles in controlling chemical weathering, while chemically weathered surface material can regulate climate change. To estimate global chemical weathering fluxes and CO₂ balance, it is important to identify the characteristics and driving factors of chemical weathering and CO₂ consumption on the Tibetan Plateau, especially in glaciated catchments. The analysis of the hydro-geochemical data indicated that silicate weathering in this area was inhibited by low temperatures, while carbonate weathering was promoted by the abundant clastic rocks with fresh surfaces produced by glacial action. Carbonate weathering dominated the riverine solute generation (with a contribution of 58%, 51%, and 43% at the QiangYong Glacier (QYG), the WengGuo Hydrological Station (WGHS), and the lake estuary (LE), respectively). The oxidation of pyrite contributed to 35%, 42%, and 30% of the riverine solutes, while silicate weathering contributed to 5%, 6%, and 26% of the riverine solutes at the QYG, WGHS, and LE, respectively. The alluvial deposit of easily weathering fine silicate minerals, the higher air temperature, plant density, and soil thickness at the downstream LE in comparison to upstream and midstream may lead to longer contact time between pore water and mineral materials, thus enhancing the silicate weathering. Because of the involvement of sulfuric acid produced by the oxidation of pyrite, carbonate weathering in the upstream and midstream did not consume atmospheric CO₂, resulting in the high rate of carbonate weathering (73.9 and 75.6 t km⁻² yr⁻¹, respectively, in maximum) and potential net release of CO₂ (with an upper constraint of 35.6 and 35.2 t km⁻² yr⁻¹, respectively) at the QYG and WGHS. The above results indicate the potential of the glaciated area of the Tibetan Plateau with pyrite deposits being a substantial natural carbon source, which deserves further investigation.  相似文献   

Numerical analysis of the effect of natural microcracks on the supercritical CO2 fracturing crack network of shale rock based on bonded particle models          下载免费PDF全文

Peihuo Peng  Yang Ju  Yongliang Wang  Shaoqin Wang  Feng Gao 《国际地质力学数值与分析法杂志》2017,41(18):1992-2013

The problem of predicting the geometric structure of induced fractures is highly complex and significant in the fracturing stimulation of rock reservoirs. In the traditional continuous fracturing models, the mechanical properties of reservoir rock are input as macroscopic quantities. These models neglect the microcracks and discontinuous characteristics of rock, which are important factors influencing the geometric structure of the induced fractures. In this paper, we simulate supercritical CO₂ fracturing based on the bonded particle model to investigate the effect of original natural microcracks on the induced‐fracture network distribution. The microcracks are simulated explicitly as broken bonds that form and coalesce into macroscopic fractures in the supercritical CO₂ fracturing process. A calculation method for the distribution uniformity index (DUI) is proposed. The influence of the total number and DUI of initial microcracks on the mechanical properties of the rock sample is studied. The DUI of the induced fractures of supercritical CO₂ fracturing and hydraulic fracturing for different DUIs of initial microcracks are compared, holding other conditions constant. The sensitivity of the DUI of the induced fractures to that of initial natural microcracks under different horizontal stress ratios is also probed. The numerical results indicate that the distribution of induced fractures of supercritical CO₂ fracturing is more uniform than that of common hydraulic fracturing when the horizontal stress ratio is small.  相似文献   

Validating the ability of the discontinuous deformation analysis method to model normal P‐wave propagation across rock fractures          下载免费PDF全文

Wenjun Gong  Yunsheng Wang  Youjun Ning  Yonghong Luo  Yujie Wang 《国际地质力学数值与分析法杂志》2017,41(10):1267-1282

The effects of fractures on wave propagation problems are increasingly abstracting the attention of scholars and engineers in rock engineering field. This study aims to fully validate the ability of discontinuous deformation analysis (DDA) to model normal P‐wave propagation across rock fractures. The effects of a single fracture and multiple parallel fractures are all tested. The results indicate that DDA can accurately reflect the fracture effects, including the fractures stiffness, the fracture spacing and the fracture number, and the effects of incident wave frequency on one‐dimensional P‐wave propagation problems. Thus, DDA is able to deal well with normal incident P‐wave propagation problems. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

Record of enhanced late Pliocene tectono‐volcanic activity in the Lesser Antilles forearc basin: Insights into the long‐term behaviour of the Lesser Antilles subduction zone          下载免费PDF全文

Jean‐Len Léticée  Mélody Philippon  Jean‐Jacques Cornée  Philippe Münch  Lyvane De Min  Mélanie Noury  Auran Randrianasolo  Jean‐Frédéric Lebrun 《地学学报》2017,29(4):245-252

In this study, we investigate the Lesser Antilles forearc basin, focusing on the late Pliocene to Pleistocene sedimentary archives in order to track the occurrence of extreme events triggered by enhanced subduction‐related tectono‐volcanic activity. We identify late Piacenzian deposits covering a major regional erosional surface, displaying sedimentary dykes and large marine boulders embedded in a mixed continental–marine matrix, characteristic of tsunamites. We interpret this episode of platform emersion and the successive cataclysmic deposits as resulting from enhanced tectonic activity at the interface of the subduction zone, synchronous with the initiation of the Lesser Antilles volcanic arc. We then discuss the implications in terms of the mechanical behaviour of the Lesser Antilles subduction zone.  相似文献   

Consolidation in spatially random unsaturated soils based on coupled flow‐deformation simulation          下载免费PDF全文

Y. Cheng  L.L. Zhang  J.H. Li  L.M. Zhang  J.H. Wang  D.Y. Wang 《国际地质力学数值与分析法杂志》2017,41(5):682-706

This paper integrates random field simulation of soil spatial variability with numerical modeling of coupled flow and deformation to investigate consolidation in spatially random unsaturated soil. The spatial variability of soil properties is simulated using the covariance matrix decomposition method. The random soil properties are imported into an interactive multiphysics software COMSOL to solve the governing partial differential equations. The effects of the spatial variability of Young's modulus and saturated permeability together with unsaturated hydraulic parameters on the dissipation of excess pore water pressure and settlement are investigated using an example of consolidation in a saturated‐unsaturated soil column because of loading. It is found that the surface settlement and the pore water pressure profile during the process of consolidation are significantly affected by the spatially varying Young's modulus. The mean value of the settlement of the spatially random soil is more than 100% greater than that of the deterministic case, and the surface settlement is subject to large uncertainty, which implies that consolidation settlement is difficult to predict accurately based on the conventional deterministic approach. The uncertainty of the settlement increases with the scale of fluctuation because of the averaging effect of spatial variability. The effects of spatial variability of saturated permeability k_sat and air entry parameters are much less significant than that of elastic modulus. The spatial variability of air entry value parameters affects the uncertainties of settlement and excess pore pressure mostly in the unsaturated zone. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

Dynamic response of a pile considering the interaction of pile variable cross section with the surrounding layered soil          下载免费PDF全文

Liu Gao  Kuihua Wang  Si Xiao  Wenbing Wu 《国际地质力学数值与分析法杂志》2017,41(9):1196-1214

Assuming that the pile variable cross section interacts with the surrounding soil in the same way as the pile toe does with the bearing stratus, the interaction of pile variable cross section with the surrounding soil is represented by a Voigt model, which consists of a spring and a damper connected in parallel, and the spring constant and damper coefficient are derived. Thus, a more rigid pile–soil interaction model is proposed. The surrounding soil layers are modeled as axisymmetric continuum in which its vertical displacements are taken into account and the pile is assumed to be a Rayleigh–Love rod with material damping. Allowing for soil properties and pile defects, the pile–soil system is divided into several layers. By means of Laplace transform, the governing equations of soil layers are solved in frequency domain, and a new relationship linking the impedance functions at the variable‐section interface between the adjacent pile segments is derived using a Heaviside step function, which is called amended impedance function transfer method. On this basis, the impedance function at pile top is derived by amended impedance function transfer method proposed in this paper. Then, the velocity response at pile top can be obtained by means of inverse Fourier transform and convolution theorem. The effects of pile–soil system parameters are studied, and some conclusions are proposed. Then, an engineering example is given to confirm the rationality of the solution proposed in this paper. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献